function y = f(x)
load('Q.mat')
y = 0.5*x'*Q*x;